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Abstract 

Effect of weak disorder on tunneling through a potential barrier is studied 
analytically. A diagrammatic approach based on the specific behavior of sub- 
barrier wave functions is developed. The problem is shown to be equivalent 
to that of tunneling through rectangular barriers with Gaussian distributed 
heights. The distribution function for the transmission coefficient T is de- 
rived, and statistical moments (T n ) are calculated. The surprising result is 

that in average disorder increases both tunneling conductance and resistance. 
PACS numbers: 73.40.Gk, 73.50.Bk, 42.25.Bs 
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It is well known that in one-dimensional disordered system the interference of propagating 
waves results in strong localization. Transmission coefficient (transmittance) T for such a 
system is not self-averaging quantity and obeys the log-normal distribution law. 

The question arises what is the effect of multiply scattering of subbarrier (evanescent) 
waves on the tunneling transmissivity, when the standard conceptions of phases, interference 
and coherent trajectories are unapplicable. This problem apparently arises in a wide variety 
of applications in solid state theory, optics and radiophysics. 

The particular question this paper is addressed to is how the one-dimensional quan- 
tum tunneling is disturbed by random perturbations of the shape of potential barrier. We 
consider a disordered region < x < L with the potential given by 

U(x) = V + v(x) , (1) 

v(x) being a random function of x. An incident (from the left) particle has energy k 2 (we 
use units where h 2 /2m = 1) that is less then the height of the unperturbed barrier Vo: 

V = k 2 + k 2 > k 2 . (2) 

Transmittance of such a system is expressed through the Green function of the problem as 

T=\r\ 2 = 4k 2 \g(0,L)\ 2 , (3) 

where r is the amplitude transmission coefficient, and Q(x,x') satisfies the equation 

(k 2 - H)g(x,x') =5{x-x'), (4) 

with the Hamiltonian 

H = ~ + U(x) [9(x) - 9(x - L)\ . (5) 

In the case Vq < k 2 there exist standard approaches based on the averaging over rapidly 
oscillating phases |TJ . In the case of subbarrier penetration evanescent waves have no rapidly 
oscillating phases, and we need another method allowing to utilize their special features. 



Below we develope a diagrammatic perturbation approach that takes an advantage of the 
rapid decay of evanescent waves. 

In order to take into account explicitly scattering by the barrier edges we decompose H 
into 3 terms so that 

H = (-^ + V ^j + AVx(x) + AV 2 (x) , (6) 

where 

AVi(x) = v(x) [9(x) - 9(x - L)] , 

AV 2 {x) = V [B(x) - 6{x — L) — 1] . 

We introduce the "subbarrier" Green function G(x,x') satisfying the equation 

/ d 2 



-AV 1 (x))G(x,x') = 5(x~x'). (7) 



\dx 2 

Then the scattering by edges described by the term AV 2 (x) can be taken into account 
through the equation 

/CO 
dx\G{x, xi)£sV2{xi)G{xi, x 1 ) . (8) 
-co 

Since AT4(a^i) = for < x\ < L no scattering take place in the region of integration in 
Eq.(H), and the behavior of both G(x, x±) and G(xi, x') is completely defined by their values 
at the edges. For x, x' G (0, L) we obtain from Eq.@ 

G(x, x') = G(x, x') - k + (G(0, x)g{0, x') + G(x, L)Q{x\ L)) , (9) 

where k± — K ± ik. It follows that 

e(0 - L) = m m-S(P(o,L)' (10) 

where 

f(z) = l + K + G(z,z). (11) 
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To proceed let us present Eq.(^) in the integral form: 

Gtx, x') — Gq(x, x') + / dxiG (x, xi)v(xi)G(xi, x'). (12) 

Jo 

Here we have introduced the unperturbed "subbarrier" Green function 

G (x, x') = — — exp (— k \x — x'\) . (13) 
2k 

The idea of further calculations can be illustrated by consideration of the first 
term G^(x,x') in the diagrammatic expansion of the average subbarrier Green function 

(G(x,x')y. 

G^(x,x') = W [ L dx 1 G (x,x 1 )G (x 1 ,x'). (14) 
Jo 

Here the correlation function of the random part of potential v (x) is taken to be 

(v(x)v(x')} = W5(x - x') (15) 
Straightforward calculation shows that 

G^\x,x') = G (x,x') + , (16) 

where 

j(x,x') = 2 - e - 2fiX - e- 2K{L ~ x '\ 

and we have introduced the subbarrier scattering length 

/ = Ak 2 /W. (17) 

The first term in brackets in Eq. (p!6|) results from the integration over the region x < X\ < x', 
where Gq(x, Xi)Gq(xi, x') does not depend on x±, while the second one comes from < xi < 
x and x' < xi < L, where Go(x, xi)Gq(xi, x') decreases exponentially as a function of 
x\. In the case that x — x' ~ L the contribution from these exponential tails is small in 
parameter (kL)~ 1 . Therefore, if kL ^> 1 the region x < x\ < x' produces the most divergent 
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with L contribution, and integration region in Eq. (|14]) can be restricted by (x, x') . This 
simple consideration being extended for arbitrary order in perturbation parameter W leads 
to conclusion that only "x-ordered" diagrams should be taken into account. This means that 
we have to allow for diagrams in which the scattering points Xi (variables of integration) 
are ordered along a straight line (Green function line) with respect to each other and to the 
fixed points x = and x' = L. For m-th order of the average Green function an ordered 
diagram corresponds to the integration over the regions 

< xi < x 2 < ■ ■■ < x m < L (18) 

if scattering points Xi are numerated subsequently from the left to the right (see Fig. la). It 
is easy to see that above arguments being extended for arbitrary order of arbitrary moments 
of Green function bring us to the following selection rules: i) only diagrams containing no 
crossings of scattering lines are important; ii) scattering points belonging to the same Green 
function line must be ordered with respect to each other. 

Selection rules being applied to the calculation of the average subbarrier Green function 
give us 

1 00 W f 

(G(x,x')) = -— exp(-n\x -x'\) ^ — / dx 1 ...dx m = (19) 

ZK m — 1 ^ K J 

m - 1 X 1 <X 2 ...<X m 

00 W L m 

= G (x,x') V — - — r = G (x,x') exp (\x - x'l //) 

m=l 4k ml 

Function (G(x, x')) contains factor describing an exponential enhancement in comparison 
with Go- It can be shown that this is also true for all higher statistical moments (G n (x, x')) 
as well. In the case kI ^> 1, kL ^> 1 arguments similar to above lead us to the necessity to 
summarize only ordered diagrams of the type shown at Fig. lb. To perform this calculation, 
let us consider a typical pair of diagrams different from each other by permutation of two 
scattering lines (Fig. 2). The sum of these diagrams can be written as 



J dy J dz... + ... J dy J dz... = ... J dy J dz. 



where we have used independence of Green functions on scattering points for ordered dia- 
grams. As one can see, the limits of integration over z do not depend on y. Repeating this 
procedure we can add other diagrams different by the location of one selected scattering 
line. This trick can be continued till we meet scattering line of the same kind (i.e., it has 
scattering points lying on the same Green function lines as that we watch for). As the result 
the integrand does not depend on coordinates, and we should no longer care of the order- 
ing between different scattering lines. Therefore, the problem is reduced to combinatorial 
governed by the fact that integration should be restricted by the scattering points of the 
nearest scattering lines. Upon applying this procedure to all possible pairs of this kind in 
all perturbation orders we deduce that for x, x' G (0, L) 



{G n (x,x')) 
(G(x,x')) n 



(G 2 (x,x')} 



c\ 



2 



(20) 



(G(x,x')} 2 

where C 2 = ^n(n — 1) is the number of possibilities to combine n Green functions into pairs. 
Calculation of (G 2 (x,x')) is straightforward and yields 

(G 2 (x, x')) = (G(x, x')) 2 exp (\x-x'\ /I) . 

It is convenient to define the quantity £ through the equation 

exp(£t/2) = G(0,L)/G (0,L) (21) 

where the dimensionless length t = L/l. From Eq. (|20|) follows 



= exp Q n ( n + ^ _ 

Statistical moments of this form correspond to the Gaussian distribution low for the random 
quantity £(£): 




ke-i) 2 



8vr 

With the same accuracy as in derivation of Eq.(|20D one can find that 

(G n (0,L)G m (x,x)) w (G n (0,L))G™(x,x). 



(22) 



This means that G(0, 0) and G(L, L) can be replaced by their unperturbed values which are 
equal to (— 1/2/c) with the accuracy 1/kI. After this procedure we obtain 

1 - cos 46* 



T(0 



(23) 



cosh(?7 — £)t — cos 46 ' 

where rj = 2kI 3> 1, 6 = tan ^(/t/fc). One can see that T(£ = 0) is exactly the transmission 
coefficient T (t) of the ideal barrier of dimensionless length t and height rj. 

1 - cos 46 



cosh rjt — cos 46 



(24) 



Eqs.(|22|) and (|23|) bring us to the conclusion that our problem is equivalent to the problem of 
transmittance of rectangular barrier with random height where deviation £ from unperturbed 
value rj is described by the Gaussian distribution function ([221). 
From Eqs.(g3),(H) follows tnat in tlie nmit t = L/l~> 1 



,L/l 



'T(L)' T (L) ' 

which coincides with the result of Ref . [Q] . 

Eqs.(|22|),(|23|) allow us to calculate all statistical moments of T 



(25) 



(T n ) 



d£P t (0T n (0- 



(26) 



For the long distance limit (t ^> 1) we obtain 



(T n ) = T n (t)e 2n ( n+ ^> , n < 77/4, t > 1. 



(27) 



The main contribution to the large statistical moments n ^> rj/4 comes from £ of the 
order of rj, corresponding to the strong fluctuations of random potential when they become 
comparable with height of unperturbed barrier. In this case our approach is not applicable 
and methods of Refs. should be used. Neglecting these rear downshooting we are able to 
retrieve the distribution function 

2N 



Pt(a 



_3 
CY * 

~ exp 



'tort 



In a 



rjt 3> In a. 



(28) 



for normalized transmittance a(t) = T/Tq (0 < a < 1/Tq) that characterizes the enhance- 
ment of the transmission rate in comparison with that of the unperturbed barrier. 

In Ref. 0] the subbarrier transmission problem was attacked with the aid of the invari- 
ant imbedding method. Due to an unestimated approximation which enabled to solve the 
invariant imbedding equations there was obtained an expression like Eq. of the present 



paper but different by term —1 in the brackets in exponent (|22|) . This difference becomes 
crucial for large t and leads to wrong distribution function of transmission coefficient. 

It should be noticed that since subbarrier scattering does not randomize phase (this fol- 
lows from the fact that all functions in Eq . fll2|) are real), our method is applicable also to cal- 
culation of the statistical moments of the amplitude transmission coefficient r = 2kQ(0,L). 
With the use of Eq. (p2|) we obtain: 

(r n ) « r n (t)e^ n(n+1) * , for n < r//2, 2kL > 1 (29) 

where To is unperturbed value. It is readily seen that even mean field characteristic (r) con- 
tains exponential enhancement factor. This fact can be easily understood by straightforward 
consideration of the Dyson equation for the true average Green function (Q(x,x')): 

(G(x,x')) = Go{x,x') + [ dxiGo(x,xi)E(xi)(G(xi,x')) , (30) 

Jo 

which can be rewritten in differential form: 

(J^ - k 2 - E(x)) (G(x u x')) = S(x - x') . (31) 

For weak disorder we can use the simplest approximation for the self-energy S(x) ~ 
WGo(x,x) which on a distance 1/k <C L becomes equal to negative constant — W/2k. 
Comparing Eq. (pTj) with Eqs. ([|), (||) one can see that disorder causes an effective lower- 
ing of the barrier height equal to its mean square deviation —W/2k. For weak scattering 
(k 2 ^> S that is equivalent to > 1) the solution of Eq. ( p0|) or ( |3lD turns out to be 
exactly the same as that obtained above by summarizing of ordered diagrams. 

To conclude, diagrammatic approach for calculation of any moments of subbarrier Green 
function has been developed. The distribution function of transmission coefficient through 
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ID disordered potential barrier has been found. It is shown that an ensemble of disordered 
barriers is equivalent to ensemble of rectangular barriers with random Gaussian distributed 
heights. Disorder causes enhancement on average of both tunneling conductance and resis- 
tance. 
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